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Abstract. 



Systems of two ordinary and partial differential equations (ODEs and PDEs) had 
been obtained from a scalar complex ODE by splitting it into its real and imaginary 
parts. The procedure was also carried out to obtain a four dimensional system by 
splitting a complex system of two ODEs into its real and imaginary parts. Systems 
of three ODEs had not been accessible by these methods. In this paper the complex 
splitting is used iteratively to obtain three and four dimensional systems of ODEs 
and four dimensional systems of PDEs for four functions of two and four variables. 
The new systems of four ODEs are distinct from the class obtained by the single 
split of a two dimensional system. Illustrative examples are provided. 



1 Introduction 

Lie [1] developed his study of the symmetry of differential equations for complex 
functions. Of course, to be different iable they have to be analytic. Though this was 
assumed, analyticity was not used in the sense of treating the complex function as 
the real and imaginary parts connected to each other by the Cauchy-Riemann (CR) 
conditions. Thus, while the complex nature of the function was very important for 
the topological properties of the Lie groups that arose, it was not used directly for the 
differential equations. More recently [2] the splitting of a scalar ODE, into its real 
and imaginary parts, was exploited to obtain methods to deal with systems of ODEs 
and PDEs. This was called complex symmetry analysis (CSA). Of particular interest 
was its application to dealing with the variational principle for systems of ODEs [3, 4] 
and to linearization (conversion of the equation to linear form by transformation of 
the dependent and independent variables) of ODEs [5, 6, 7]. The latter was of special 
interest because of the connection between geometry and the symmetries of systems 
of ODEs [8, 9, 10, 11]. This enabled one to use geometric methods to linearize 



systems (including the scalar case) of ODEs [12, 13]. This method allowed one to 
not only write down the linearizing transformation but also to directly provide the 
solution. The procedure also led to some new insights regarding standard systems 
of equations and to methods for solving systems that were not amenable to solution 
by standard symmetry analysis [4, 14]. 

One might have expected that a system of four ODEs or PDEs could be obtained 
by using quaternions. This turns out to be impossible. The reason is that while 
algebra works for quaternions calculus does not. Consider q = w + ix + jy + kz, 
subject to the usual quaternion rules that i 2 = j 2 = k 2 = — 1 , ij = k = —ji ,jk = 
i = —kj ,ki — j — —ik. For dq/dq = 1 the derivative operator is 

d d d . d , , d /1X 
— = a— + 10^- + jc— + kd— , (1) 
dq aw ox Oy Oz 

subject to the condition that a — b — c — d = 1. Now ask that dq 2 /dq = 2q. This 
requires the above condition along with the additional conditions, a — b = a — c = 
a — d = 1, which are inconsistent with the earlier requirement. One has to obtain 
the four dimensional system by other means, such as using a complex system of two 
ODEs and splitting it into a system of four ODEs. 

Though one can split a system of two complex ODEs into one of four ODEs by this 
procedure, it is not possible to use it to obtain a system of three ODEs from either 
a scalar or a vector equation. In fact, all odd dimensional systems are inaccessible 
by the usual CSA splitting procedure. However, one would like to be able to use the 
power of CSA for three (and other odd dimensional systems) dimensional systems 
as well. In this paper the complex splitting is used iteratively to be able to generate 
three and four dimensional systems of ODEs and to obtain systems of four PDEs 
for four functions of two or four variables. The procedure could be used with more 
than two iterations to generate higher dimensional systems as well, but we will not 
follow that up here. 

The method used to obtain the system of ODEs is to start with a scalar ODE and 
regard the dependent variable as a complex function of a real variable, as was done 
for CSA, and split it into a system of two ODEs. Now regard the two dependent 
variables as themselves complex functions of a real variable. This provides a system 
of four ODEs. Of course, if both steps had been merged into one, only a system of 
two equations would again have been obtained. (This would be the reduced system 
obtainable from the new system of four equations.) Instead, we first close our eyes 
to the fact that we are going to treat the two dependent variables as complex, and 
only after obtaining the system of two equations do we treat each of the dependent 
variables as themselves complex. The result of repeating the CSA procedure is 
different from treating the dependent variable of the original scalar equation as a 
function of two complex functions of a real variable as the symmetry structure of 
the systems is different. 
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In the CSA procedure it was found that one may be able to linearize the base 
scalar equation even though the corresponding real system does not have enough 
symmetries to allow linearization, or even to permit of solution by usual symmetry 
methods. The system of four ODEs could then have much fewer symmetries while 
the base system is linearizable, or even have a base system with too few symmetries 
while the scalar equation is linearizable [14]. It is to be expected that here we will 
not only get the same sort of cases, but that there would be some even stranger 
cases arising than for CSA. 

One has another option. For the system obtained by the first split, one can treat 
one of the two dependent variables as real and the other as complex. This provides 
the system of three ODEs. In fact, it provides two systems of three ODEs, as we can 
choose either of the dependent variables to be real and the other complex. The two 
systems obtained are dual to each other in some sense. That sense will be clarified 
by some examples that we will provide. 

Of course, one could take the CSA system of PDEs and treat each of the dependent 
variables as complex functions of the two real variables. This would yield a system of 
four PDEs for four functions of two variables. Alternatively, we could have started 
with the system of two ODEs and now treated both the dependent and independent 
variables as complex. This split also gives a system of four PDEs for four functions 
of two variables. The two systems obtained are dual to each other in a fairly obvious 
way. Splitting the system of two PDEs by treating the dependent and independent 
variables as complex gives a system of four PDEs for four functions of four variables. 

The role of the CR-equations was fairly obvious in the original CSA and had not 
been spelled out explicitly. A geometric description of these equations was given 
later [14] but the detailed requirements for a general scalar ODE were not. On 
the double use of the splitting procedure the conditions become more thoroughly 
coupled and need to be spelled out explicitly. It turns out that even for the original 
CSA the conditions for the derivatives of the functions involved on the right side of 
the semilinear equations with respect to the derivatives of the dependent variables 
are more complicated than was envisaged. These have been stated explicitly here. 

We have limited the discussion to second order equations only. One can, of course, 
go to higher order equations but that complicates the expressions without providing 
any further understanding of the procedures being developed. Also, there is no 
direct equivalent of the geometric connection and procedures for the higher order 
equations. (There is an indirect connection provided by differentiating the second 
order ODE and requiring that the original equation hold [15, 16], but that will not 
be followed up here.) On the other hand, one could have limited oneself to first 
order equations but that is the degenerate case and results for it will not hold more 
generally. 
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The basic view adopted here is to use the connection of a second order system 
of differential equations with a scalar second order ODE to be able to solve, or 
otherwise deal with, the given system. For this purpose, one ideally needs clear 
criteria to be able to determine whether the given system is, or is not, related to 
some scalar ODE. The criteria should be such that a computer code could be written 
to check the given system for the relationship, and if it is related to construct the 
required ODE, which could then be appropriately dealt with. This had not been 
done even for the original CSA. Here we state the explicit criteria as theorems for 
the original CSA and then for the double split systems. 

The plan of the paper is as follows. In the next section we present the basics of the 
CSA splitting procedure and briefly mention symmetries of differential equations. 
We also present two characterization theorems there that had not been provided 
earlier. In the subsequent section we give the split into a system of three ODEs. In 
section four we give the split into a system of four ODEs. The section after that 
deals with the split into a system of four PDEs for four functions of two variables. 
In section six we present the split into a system of four PDEs for four functions 
of four variables. For each of the systems we provide some examples in the same 
sections. In the concluding section we give a brief summary and discussion of the 
results. 



2 Complex Splitting and Symmetries 

Consider a general second order ODE 



We can now take either u to be a complex function of the real variable r or take 
both u and r to be complex. Let u = p + iq in the former case and put 



u"(r) = f(r;u,u') . 



(2) 



f(r; u, u') = f r (r; p, q; p', q') + if l (r; p, q; p', q') . 



(3) 



The resulting system of ODEs is 



P"{r) = f(r;p,q;p',q') ,q"(r) = f(r;p,q;p',q') . 



(4) 



The CR-equations for this system are 



p J q i J q J p 1 J p' J q' 1 J q' J p' 1 



(5) 



with no conditions on p and q other than second differentiability. 



For the latter case let r = s + d as well. Then 



f(r;u,u') = f r (s,t;p,q;p s ,q s ,pt,q t ) + tf(s,t;p,q;p s ,q s ,p t ,q t ) . 



(6) 
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Splitting the equation into its real and imaginary parts then gives the system of two 
PDEs 

Pss - Ptt + 2q st = 4/ r , q ss - q tt - 2p st = 4f . (7) 

The CR-equations for this system include conditions for p and q, apart from the 
previous ones (which are now more complicated), 

Ps = Qt , Pt = -q 8 ; 

/r £% rr £% . 

S Jt 1 Jt J S 5 

/; = r q , r q = -r P ■ (») 

The derivative of the functions with respect to the derivatives is more complicated 
to state. The problem is that 

Thus, for /(r; u, u') to be analytic, f r and f l cannot depend arbitrarily on p s , q s , p t 
and qt but must depend on p s + q t and pt — q s . If we call these variables and ip, 
respectively, then the last conditions are 

fs = fibi fib = ~f<b ■ (10) 



Though the CR-equations were taken as obvious, they are needed to characterize 
systems of (real) ODEs and PDEs that can arise by splitting a scalar ODE by CSA 
methods. To complete the CSA procedure we state the following two theorems. 



Theorem 1: A system of two second order ODEs (4) corresponds to a scalar second 
order ODE (2) if and only if it satisfies the system of CR-conditions (5). 



Theorem 2: A system of two second order PDEs (7) corresponds to a scalar second 
order ODE (2) if and only if it satisfies the system of CR-conditions (8) and (10) 
where (f> = p s + q t and ip = p t — q s . 



Symmetry analysis deals with the infinitesimal generators that leave the differential 
equation invariant under point transformations (r, u) — > (R, U) say, 

X = ^ r ,u)^ + V (r,u)^-. (11) 

To be able to apply the operators to differential equations these generators have to 
be prolonged or extended to include the higher derivatives. For the first derivative 

x,11 =4+'l;+' [11 ^"'^ < 12 > 

where 
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^[o] _ ^ ^i^ r s t anc i s f or the total derivative, 



d_ 

dr 




9 , 9 , „." 9 



(14) 



<9r du du' 



and -u n (r) stands for the n th derivative. The prolonged generators for higher order 
differential equations can be similarly obtained using (13). 

We can extend the analysis to systems of ODEs by replacing the scalar u by a vector 
u and the corresponding partial derivative by V u . Consequently we must replace 
the scalar r\ and rj 1 ^ by the vectors r\ and The extension to the PDEs is more 
complicated. The scalar r has to now also be replaced by a vector s and along with 
it the scalar £ by the vector £, but now the derivative of u becomes V s u and the 
derivative with respect to this vector of partial derivatives becomes too messy to 
read. As such, we write V s u = ui. Of course, we need to also bear in mind that the 
variables that the functions depend on will not be the derivatives given but linear 
combinations as we saw for the CR-conditions. The second derivative can then be 
written as u 2 and so on. For the complex case the X was written as Z and similar 
changes were made for the coefficients but it will be more convenient to use the same 
notation throughout here. 



3 System of Three ODEs by Double Splitting 



Consider (2) as the base scalar equation with the split of the function given by (3) 
leading to the system (4) subject to the CR-equations (5). Now regard p as a real 
variable x and q as the complex variable y + iz. We run into a problem here. There 
are three second order ODEs but the number of functions to be obtained from / 
must be even. To avoid this problem we put 



The choice of what part to put into g and what part to put into G is clearly arbitrary. 
For definiteness, we define g to consist of all those terms that do not involve q or q 1 . 
Then G consists of all those terms that do. Now we proceed with the second split 
by putting 



f(r; u, u') = g(r; p, p') + t G(r; p, q; p, q) . 



(15) 



G(r; p, q; p, q') = k(r; x, y, z; x', y', z') + i l(r; x, y, z; x', y', z') , 



(16) 



yielding the system of three ODEs 



x" = h(r; x, x) , y" = k(r; x, x') , z 



/(r;x,x) , 



(17) 



subject to the CR-equations 




(18) 
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The first prolonged symmetry generator is: 

+rf W(r; x, V x ,)^ + ^(r; x, V*)^ + rf^r; x, V x ,)^ . (19) 



Instead of this procedure, at the second split we could have taken p(r) = x(r)+i y{r) 
and set q(r) = z(r). This will give a dual system in some sense. We would now have 
to set 

f(r; u, v!) = g(r; x, y, z; x\ y', z') + i h(r; x, y, z; rf, y', z') + t k(r; z, z') , (20) 

and the slightly modified CR-equations 

9y = h z , g z = —h y ; g y i = h z > , g z > = —h y i . (21) 

The prolonged symmetry generator remains unaltered in form except that now the 
coefficient of the last term is pure imaginary, to account for the last imaginary term 
in (20). The sense of the duality will be clarified by the examples. 

For completeness we state a theorem for the characterization of systems of three 
ODEs that correspond to a scalar ODE by double splitting. 

Theorem 3: The system of three ODEs (17) corresponds to the scalar ODE (2), 
for any consistent identification of the function given by (15) and (16), provided the 
CR-conditions (18) hold. 

Remark: The dual procedure gives the same system but now we require that the 
split given by (20) and (21) holds. 

Example 1: Consider the free-particle scalar ODE, u" = 0. It clearly yields the 
system x" = y" = z" = 0. The splitting of the functions is obviously trivial. How- 
ever, the infinitesimal symmetry generators of the system are not trivially related to 
the generators of the scalar ODE. Even for the original CSA, it had been noted that 
the symmetries for the system could not be a simple doubling of the symmetries 
of the original ODE, as the maximal Lie algebra for the system is s/(4,R), which 
has 15 generators, while the algebra for the scalar free particle equation is sl(3, 1R ), 
which has 8. Doubling gives one extra generator. It is also clear that it cannot be 
a simple matter of leaving one generator out, as the complex generators occur in 
pairs. What happened there was that we lost two generators and got one new one. 
The system must clearly have 24 generators, as the algebra is s/(5,!R). However, 
the double splitting gives only 23 operators of which 15 are symmetries and 8 are 
Lie-like [2], five of the nine dilations coming from the dependent variable and the 
four local projective symmetries are missing. We have to use the actual symmetries 
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obtained here and require closure of the algebra to generate the full 24. The dual 
system also gives the same symmetry structure, as expected. However, due to the 
fact that the k in (20) has a coefficient with iota, even though k itself is zero, the 
operator for the symmetry carries an imaginary. If we do not put in the iota for the 
operator, we lose some of the symmetry generators. 

Example 2: Consider the scalar equation [17] 

u" = s" V . (22) 

It splits into the following system of ODEs: 

x" = -2s~ 5 yz ; 
y " = -2s~ 5 zx ; 

z " = -2s- 5 xy ; (23) 

subject to the further algebraic constraint 

X 2 + y 2 = z 2 . (24) 

There are no symmetries among the 8 Lie-like operators obtained from the splitting 
of the complex generators 

d d d d 

Xi = s— + 3u— ; X 2 = s 2 — + su— . (25) 

OS OU OS ou 

However, the system admits the two symmetry generators: 

v d d d 

ox oy Oz 

Y 2 = s 2 — + sx— + sy— + sz— . 26 

os Ox Oy Oz 

Example 3: Consider the scalar (Emden-Fowler) equation [18] 

u" + 5s _1 «' + u 2 = . (27) 

It splits into the following system of ODEs: 

x" + 5s - V - 2yz = ; 
y" + 5s _1 y' — 2zx = ; 

z" + 5s - V - 2xy = ; (28) 
subject to the further algebraic constraint 

X 2 + y 2 = z 2 . (29) 
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The Emden- Fowler equation has only the one scaling symmetry 

x - 4 S - 2u l ■ (3o) 

The split system has 4 Lie-like operators, none of which are symmetries of the 
system. However, the system admits the scaling 

d d d d 

Y = s-^ — 2x— - 2y— - 2z- . (31) 
os ox Oy Oz 



4 System of Four ODEs by Double Splitting 

For the four dimensional system, after the first split of (2) given by (3), (4) and (5) 
we could set p(r) = w(r) + i x(r) and q(r) = y(r) + i z(r) and 

f(r; p, q; p', q') = g(r; w, w') + i h(r; w, w') , 
fir; p, q; p', q') = k(r; w, w') + t Z(r; w, w') , (32) 

yielding the system of four ODEs 

w"(r) = g(r; w, w') , x"(r) = h(r; w, w') , 

y "(r) = k(r; w, w') z"(r) = l(r; w, w') , (33) 

subject to the CR-conditions 

9w h X ky ~\- l Z , Cj X H w k z ly , 

9y kyj , (J Z Hy k X ~\~ lyj , 

g w i + h x i = kyi + l z i , g x i h w / = k z > lyi , 
g y t + h z i = —k w i — l x i , g z i — h y i = —k x i + l w > , (34) 

where w = (w,x,y,z). The prolonged symmetry generator can now be written as 

d 

x = £( r ' w )^ + iz( r > w )- v w + !z [ %; w V)-v W ' . (35) 

Writing this equation out in detail makes it too unwieldy to convey much wisdom. 

There is no dual system to this as the splitting is direct. We could have obtained a 
system of four ODEs by a three stage splitting, setting one of the dependent variables 
in the systems of three ODEs as a complex pair. For each of the three dimensional 
systems obtained, one splitting would give the above system and one would give a 
new system. We are, here, limiting our discussion to a two-step splitting only. 

The characterization theorem is: 

Theorem 4: The system of four second order ODEs (33) corresponds to the scalar 
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ODE (2) by double complex splitting if and only if the CR-conditions (34) hold with 
the splitting (32). 



The examples will illustrate our procedure further. 

Example 4: The free particle scalar equation obviously yields the free-particle 
system of four equations, w" = x" = y" = z" = 0. The symmetry algebra must be 
s/(6, R ), which has 35 generators. There are a total of 18 symmetries and 8 Lie-like 
operators. Again, the missing ones come from the dilations involving the dependent 
variables and local projective symmetries. Again, the closure of the algebra starting 
with the derived symmetries generates the full s/(6,R). 

Example 5: Consider (22) and now split into the system of four ODEs: 

w " = s ~ 5 (w 2 - x 2 - y 2 + z 2 ) ; 
x" = s' 5 {2wx - 2yz) ; 
y " = s ~ 5 (2wy - 2xz) ; 

z" = s~ 5 (2wz + 2xy) . (36) 

This has 8 Lie-like operators of which none are symmetries of the system. However, 
the system does admit the two symmetries 

d d d d d 

Yi = s— + 3w— + 3a;— + 3y— + 3z— ; 
os aw ox ay oz 

Y 2 = s — + sw— + sx— + sy— + sz— . (37) 
os ow ox oy oz 



Example 6: Again consider the scalar (Emden- Fowler) equation (27). It splits into 
the following system of ODEs: 



w" + Ss' 1 ™' + w 2 



x 2 - y 2 + z 2 = 
x" + 5s~V + 2wx - 2yz = 
y" + hs^y' + 2wy -2xz = 
z" + hs~ l z' + 2wz + 2xy = 0. 



(38) 



The split system again has 4 Lie-like operators, none of which are symmetries of the 
system. However, the system admits the scaling 



,,9 d d d d 

Y = s- 2w— - 2x— - 2y— - 2z— 

os ow ox oy oz 



(39) 
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5 System of Four PDEs for Four Functions of Two 
Variables by Double Splitting 



For the PDEs, we need to treat the dependent variable as complex. However, there 
are two ways of doing so if we are to obtain functions of only two variables. We could 
retain the independent variable as real in the first step and then make it complex in 
the second step, or first treat it as complex and then retain the same independent 
variables in the second step. We follow the former procedure first. Starting with 
(2), with the first step splitting given by (3), (4) and (5), we put r = s + t t and 
proceed with treating the dependent variables as complex exactly as in the case for 
the system of four ODEs. 

The system of equations is now 

w ss - w tt + 2x st = 4g(s; w; w s , w t ) , x ss - x tt - 2w st = 4h(s; w; w s , w t ) ; 

y ss - Vtt + 2z st = 4fc(s; w; w s , w t ) , z ss - z tt - 2y st = 4Z(s; w; w s , w t ) , (40) 

subject to the CR-equations 

w s — x t , w t — —x s , Us — z t j Vt — ~ z s ; 
g w h x , g x h w , gy h z , g z hy , 

kw l>x i k x lw i ky ■> k z l w . (41) 

Now, as for (10), defining 

= Ws + Xt , i\) = w t - s x , K = y s + z t , X = y t - z s , (42) 
the conditions for the derivatives with respect to the derivatives can be written as 
g<f> = h^p , g^p = —htf, ; k K — l\ , k x = —l K ■ (43) 



Writing (s, t) = s, for the infinitesimal symmetry generator we will now have two 
components for ^, namely (£ s , ^*) and 7] will have four components as for the system 
of four ODEs. The additional feature is that the prolonged derivatives will be for 
V Ws and the coefficients will be Thus the prolonged generator can be written 

as 

X = e(s,w).V s + r Z (s,w).V w + r Z [11 (s,w).V Vs w , (44) 
where rj v w(s, w) is the generalization of 77M for the case of PDEs. 

The characterization theorem here is: 

Theorem 5: The system of four PDEs for four functions of two variables (40) 
corresponds to the scalar ODE (2) by double complex splitting, if and only if the 
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CR- conditions (41) - (43) hold, provided that g, h, k, I depend on the derivatives only 
in the combinations given by (43). 



Instead, if we had put r = s + t t in the first step and then proceeded to split the 
dependent variables twice, we would have got the dual system 

w ss - w tt + 2y st = Ag(s; w; w x , w t ) , x ss - x tt + 2z st = 4/i(s; w; w x , w t ) ; 
Vss - Vu - 2w s t = 4A;(s; w; w x , w t ) , z ss - z tt - 2x st = 4/(s; w; w x , w t ) . (45) 

The CR-conditions are considerably more involved. The simple ones are 



w s = y t , w t = -y s , x s = zt , x t = -z s ; 

g-w h x ky -\- l z , g x h w k z ly , 
g y + h z = —k w — l x , g z — h y = —k x + l w . (46) 

For the derivatives with respect to derivatives we have to define the new variables 

a = w s + y t , f3 = x s + z t , 7 = w t - y s , 5 = x t - z s (47) 

to get 

9a + hp = kj + l s , gp - h a = k s - Z 7 , 
9l + h s = —k a -lp , g 8 -h 1 = -kp + l a . (48) 

The form of the symmetry generator remains unchanged. Once again, we rely on 
the examples to illustrate the procedure. 



Here we need a separate theorem because the systems are apparently different, 
though they are dual to each other in some sense. 

Theorem 6: The (dual) system of four PDEs for four functions of two variables 
(45) corresponds to the scalar ODE (2) by double complex splitting, if and only if the 
CR-conditions (46) - (48) hold, provided that g, h, k, I only depend on the derivatives 
in the combinations given by (48). 



Example 7: Consider the free particle equation split into the system of four PDEs 
for two independent variables 

w ss - w tt + 2x st = ; 
x ss ~~ x tt ~~ 2w st = ; 
Vss - ytt + 2z st = ; 

z ss - z tt - 2y st = . (49) 



The symmetry generators split into 28 Lie-like operators of which 20 are symmetries. 
However, the system admits infinitely many symmetries. 
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The dual system is very similar. More precisely, it is the the system (45) with the 
right side set equal to zero. Since the CR-conditions are trivial there is no significant 
difference between the original and the dual system. 



Example 8: Again consider the scalar (Emden- Fowler) equation (27). It splits into 
the following system of four PDEs of two independent variables: 

w ss - w tt + 2x st + 10 g2 (w a + x t ) + 10 g2 1 2 (x s - w t ) 

+4{w 2 - x 2 - y 2 + z 2 ) = ; 
x ss - x tt - 2w st + 10 — ^ 2 (x s - w t ) - 10 ^ 2 t 2 (w s + x t ) + 8(wx -yz) = ; 

s t 
Vss - Vtt + 2z st + iO g2 + t2 {Vs + z t ) + io g2 + t2 ( z s " Vt) + K W V -xz) =0 ; 

z ss - z tt - 2y st + 10— ^— (z a - y t ) - 10^-^y, + z t ) + 8(wz + xy) = . (50) 

This split system again has 4 Lie-like operators, none of which are symmetries of 
the system. However, it admits the scaling 

d d d Odd 

as ot aw ox oy oz 



6 System of Four PDEs for Four Functions of 
Four Variables by Double Splitting 

This is the most straight forward (and the most complicated) of the various possibil- 
ities considered. At the first step we regard both the independent and the dependent 
variables as given by (6) to (8). For the second step we run short of symbols for the 
variables. As such, we now write the first variable (previously written the 
complex variable s + i t and the second variable (previously written as t), as the 
complex variable u + t v and write s for (s, t, u, v). Further, we put 

p(s,t) ->■ w(s) + i x(s) , q(s,t) ->■ y(s) + i z(s) , (52) 

f r (s,t;p,q;p s ,q s ,p t ,q t ) = g(s; w, V s w) + t h(s; w, V s w) ; 
f l (s, t; p, q; p s , q s , p t , q t ) = fc(s; w, V s w) + i Z(s; w, V s w) . (53) 

The system of equations is 

w ss — w tt + 2x st — w uu + w vv — 2x uv + 2y su — 2y tv 
+2z sv + 2z tv = Ag(s; w, V s w); 
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x ss x tt 2w s t %uu x vv 2w uv -|- 2z su 2Zf v 

-2y sv - 2y tv = 4/i(s; w, V s w); 
Uss ~ Utt + 2z st — y uu + y vv — 2z uv + 2w su — 2w tv 
+2x sv + 2x tv = 4fc(s; w, V s w); 

-2w sl) - 2w tv = 4/(s; w, V s w); 

subject to the CR-conditions 

w s + x t = Du + Z V ) W t ~ X S = Uv ~ Z U J 

w u -\- x v t/ s , w v x u y v -\- z u , 
9s ht = k u -\- l v , gt h s = k v — l u , 
9u ~\~ h v = —k s — It i 9v ~~ h u = —k t + Z s ; 

fi'ui "I" ky ~\- l Z , H w k z ly , 

fi 1 ?/ hz ^ w ^ x i 9z hy k x -\- l w . 
The derivatives with respect to the derivatives require the variables 

a = w s + x t + y u + z v , (5 = w t - x s + y v - z u ; 
7 = w u + x v - y s - z t , 5 = w v - x u - y t + z s . 

Then the rest of the CR-conditions are 

9a - hp = fc 7 - l s , gp + h a = k s + lj ; 
g 1 -hs = -k a + lp , g s + h^ = -hp - l a . 

The characterization theorem in this case is: 
Theorem 7: The system of four second order PDEs for four functions of four vari- 
ables (54) corresponds to the scalar second order ODE (2) by double complex splitting 
provided the functions g, h, k, I depend on the derivatives only in the combinations 
given by (57) and the CR-conditions (55) and (57) hold. 

The prolonged symmetry generator for the system is 

XW =e(s,g).V s + !Z (s,g).V g + !Z [1] (s,g,V s g).Vv s g • (58) 

The derivatives with respect to derivatives are to be taken bearing in mind the dis- 
cussion for the CR-equations. However, even if we ignore it in taking the derivatives, 
no error will ensue. 

We again rely on the examples to illustrate our systems. 

Example 9: The free-particle system of equations is given by (54), with the right 
side set equal to zero. The CR-conditions are trivial. There are now 32 Lie-like 
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(55) 
(56) 
(57) 



operators, of which only 24 are symmetry generators. As before the local projective 
symmetries are lost. Here there are 8 such. However, the dilations are not lost here. 
The system, itself, has an infinite number of symmetry generators. 

Example 10: Consider the double splitting of the Emden-Fowler equation (27) into 
a system of four PDEs for four functions of four variables. It is given by (54) with 

g = C[(sA + tB)a + (tA - sB)(5 + (uA + vB)~f + (vA - uB)5] 

-w 2 + x 2 + y 2 - z 2 ; 

h = C[(sA + tB)/3 - (tA - sB)a + (uA + vB)5 - (vA - uB)j\ 

—2wx + 2yz ; 

k = C[(sA + tB)~f + (tA - sB)5 - (uA + vB)a - (vA - uB)/3] 

—2wy + 2xz ; 

I = C[-(sA + tB)8 + (tA - sB)-i + (uA + vB)(5 - (vA - uB)a] 

-2wz - 2xy ; (59) 

where 

C = - - 5 - , A = s 2 - t 2 + u 2 - v 2 , B = 2st + 2uv . (60) 
A 2 + B 2 

There are 4 Lie-like operators none of which are symmetries. The system, like the 
scalar equation, has a scaling symmetry. 



7 Summary and Discussion 



In this paper we have considered systems of three and four second order ODEs, and 
systems of four second order PDEs for four functions of two or four variables, that 
correspond to a scalar equation, that we shall call a base equation by a specific pro- 
cedure, that we call double complex splitting. We have also provided characterization 
criteria for such systems to correspond to the base equation and a clear procedure 
to be able to construct the base equation. Thus, in principle, we could write a 
computer code that could take any such system given and check if it corresponds to 
a base equation. It could then construct the base equation. 

What is the advantage of having such base equations and constructing them? The 
point is that it is much easier to deal with the base equation than with the system. 
Thus, for example, if the base system has two infinitesimal symmetry generators 
it could be solved by symmetry methods. In fact it could have eight symmetry 
generators and thus be linearizable. In that case we could write down the solution 
directly. Following the double splitting procedure by which the system corresponds, 
we could then write down the solution for the system of ODEs or PDEs. Note that 
in this procedure the system need not have the required symmetry for being directly 
solvable by symmetry methods. 
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The procedure adopted for CSA could only give an even dimensional system as it 
simply split n (complex) equations into In (real) equations. Also, the number of 
independent variables in the PDE equals the number of dependent variables. In 
the complex double split procedure we are being considerably more adventurous. 
Having obtained the system of two "real" equations we conveniently forget that 
they arose from a scalar complex equation, treat it as complex and promptly split 
the equations again. Now we have the earlier freedom of choosing the independent 
variable to be either real or complex, while treating the dependent variables to be 
complex, but we have the additional freedom to choose to treat one of the dependent 
variables as real and the other as complex. This provides the possibility of obtaining 
an odd dimensional system. Further, to obtain the PDEs, we could choose to treat 
the independent variables as real first and then complex, complex first and then real 
or complex both times. Thus we also get the system of four PDEs for four functions 
of two variables. 

In the cases of full double complex splitting, where either the dependent variables 
were fully split or both the dependent and independent variables were double split, 
giving the system of four ODEs or four PDEs for four functions of four variables, 
there were no complications of additional dual systems arising. However, in the 
case of the system of three ODEs or the system of four PDEs for four functions of 
two variables, we got dual systems arising. The duality in the former case was very 
obvious but in the latter it was considerably more involved on account of the CR- 
conditions. Even in the case of the system of three ODEs the symmetry generators 
had to carry an iota to make sense. The "duality" of these systems needs to be 
better understood. Note that a in (56) is simply V s .w. It would be interesting to 
find out what the operators for the other variables are. Presumably, they would 
be "dual" divergence operators in some sense. This may shed some light on the 
structure of the double split systems. 

For the system of three ODEs obtained by double splitting, it would be of interest 
to consider the ambiguity due to the choice of the function g in (15), or k in (20) 
for the dual system. In some sense all choices must be "equivalent" . The question 
is whether one gets an equivalence class. Further, would they be equivalent under 
point transformations or possibly some more general transformations like contact or 
higher order transformations [19]. 

The algebraic constraint that arises in the system of three ODEs was not apparent 
in setting up the system but was found in the examples. It is interesting to note 
that it geometrically amounts to the solution lying on a cone. This seems to be 
generic for the three dimensional system. It also shows that it can be written as a 
system of two ODEs. However, that system is much more complicated. 

One would have hoped that for the system of PDEs corresponding to a base ODE 
one could use the symmetries of the ODE to obtain a "core" set of symmetries for 
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the system of PDEs. Even if the system has infinitely many symmetries, the base 
equation can only have a finite number. However, the examples show that we can 
lose all the Lie symmetry generators and be left only with Lie-like ones. This applies 
even to the free particle equation. Further for the Emden-Fowler equation, we are 
left with no Lie symmetries from the Lie-like operators, though the equation has a 
scaling symmetry and so does the double-split system. In general, we obtain Lie- 
like operators and not Lie- symmetry generators that would form an algebra. The 
Lie-like operators somehow encode the symmetries of the base equation. It would 
be most important to learn how they do so. It may be that the CR-conditions will 
enable us to re-construct the Lie from the Lie-like symmetry. 

It is of interest to note that not only for the PDEs but also for the systems of ODEs, 
we get Lie-like operators arising and lose some Lie-symmetry generators. It would 
be worth while to see this encoding of symmetry as distinct from the PDE case. 

We hope that in future the use of these systems for the variational principle and with 
linearization will be followed up. It should lead to interesting and useful results. 
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